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Abstract 

Angiogenesis is a key process in the tumoral growth which allows the cancerous tissue to impact on 
its vasculature in order to improve the nutrient's supply and the metastatic process. In this paper, we 
introduce a model for the density of metastasis which takes into account for this feature. It is a two 
dimensional structured equation with a vanishing velocity field and a source term on the boundary. 
We present here the mathematical analysis of the model, namely the well-posedness of the equation 
and the asymptotic behavior of the solutions, whose natural regularity led us to investigate some 
basic properties of the space Wdiv(fi) = {V £ L ; div(GV) 6 L 1 }, where G is the velocity field of 
the equation. 
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1 Introduction 



In the seventies, Judah Folkman puts forward the assumption that a cancer tissue, like other tissues, 
needs nutrients and oxygen conveyed by the blood vessels. Consequently, tumoral growth and develop- 
ment of metastasis are dependent on angiogenesis, a process consisting in building and developing the 
vascularization. From this discovery, a new anti-cancer therapeutic way is open : to starve cancer by 
depriving it of its vascularization. If for the last two decades, more than ten antiangiogenics drugs have 
been developed, mainly monoclonal antibodies and tyrosin kinase inhibitors, the administration proto- 
cols are far from being optimal. It is enough for example to consult the publication [13] to realize the 
paroxystic effects they can induce. 

Thus, a tool for in silico studying the administration protocols for antiangiogenic drugs could largely 
contribute to optimize the effectiveness of the treatments, in particular to avoid some therapeutic failures. 
In this direction, the construction of a mathematical model taking into account the mechanisms of tumoral 
angiogenesis and the effect of the antiangiogenics agents proves to be an essential stage in order to improve 
the use of antiangiogenic therapies. Some work (for example in [22] and [5]) was made with the aim of 
qualitatively studying the effects of antiangiogenic therapies on the control of the primitive tumor growth. 
In this work we propose a modeling which purpose is to describe the action of the currently used clinical 
protocols, not only on the tumoral growth, but also on the production of metastases. 

The model is a combination of the PDE model for the metastasis density proposed by [17] and studied 
in [3, 9], with the ODE model for each metastasis' growth of Folkman ct al. [16]. This transport equation 
endowed with a non-local boundary condition expressing creation of metastasis can be classified as part 
of the so-called structured population equations arising in mathematical biology which have the following 
general expression 

( a t p + div(F(t,x,p)) = -n(t,x,p) n 

(1) I -G-~Jp{t 1 a)=B{t,<T 1 p) a e dfl s.t. G ■ ^(a) <0 . 

{ p(0,X) = p°(X) Q 

The introduction of such equations in the linear case is due to Sharpe and Lotka in 1911 [24] and 
McKcndrick in 1926 [18]. Although these equations have been widely studied both in the linear and 
nonlinear cases (for an introduction to the linear theory see the book of Perthame [20] and to the 
nonlinear one see the book of Webb [27], as well as [21] for a survey), a complete general theory has not 
been achieved yet, even in the linear case. Indeed, most of the models have the so called structuring 
variable X being one-dimensional and often representing the age, thus evolving with F(t,a,p) = p. A 
difficulty on the regularity of solutions is introduced when the velocity is non-constant and vanishes (see 
[3, 9]). Dealing with situations in dimensions higher than one is not a common thing. 

In our case, the model is a linear equation, with 

F(t,X,p) = G(X)p 

structured in two variables : X = (x, 9) with x the size of metastasis and 9 the so-called "angiogenic 
capacity". The velocity field G vanishes on the boundary of the domain, which is a square. Moreover, we 
have an additional source term in the boundary condition of the equation : 

-G ■ l?p(t, a) = N(a) J P(X)p(t, X) + f(t, a). 

As far as we now, the mathematical analysis for multi-dimensional models is done only in situations where 
one of the structured variables is the age and thus with the first component of G being constant (see 
for instance [26, 1, 11]). In the context of the follicular control during the ovarian process, a nonlinear 
model structured in dimension two with both components of the velocity field G being non-constant is 
introduced in [14] but no mathematical analysis is performed due to the complexity of the model. 
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In the present paper, we address the problem of the mathematical analysis of our model, namely : 
existence, uniqueness, regularity and asymptotic behavior of the solutions. Following the method used in 
[2] and [3], we use a semigroup approach to deal with the existence and regularity of the solutions. The 
main difficulties we have to deal with in this two dimensional problem come from the singularity of the 
velocity field, as well as the presence of a time-dependent source term in the boundary condition. During 
the study, we take a particular attention on the problems of regularity of the solutions and approximation 
of weak solutions by regular ones, which led us to study the space W<iiv(ft) (see the appendix) . The paper 
is organized as follows : in the section 2 we present the model, in the section 3 we study the properties 
of the underlying operator and in the section 4 we apply our study to the evolution equation from our 
model. 



2 Model 

The model we developed is an improvement of the model proposed by [17] and studied in [3]. We want 
now to take into account the key process of angiogenesis in the tumoral growth and integrate it in the 
metastatic evolution. To do this, we combine a renewal equation describing the evolution of the density 
of metastasis with an ODE model of tumoral growth including angiogenesis developed by Hahnfeldt et 
al. in [16]. 



2.1 The ODE model of tumoral growth under angiogenic control 

We present now the model of Hahnfeldt et al. from [16]. Let x(t) denote the size of a given tumor at 
time t. The growth of the tumor is modeled by a gompcrtzian growth rate, which expression is : 



(2) gi{x) = axln 



where a is a parameter representing the velocity of the growth and 9 the carrying capacity of the envi- 
ronment. The idea is now to take 9 as a variable of the time, representing the degree of vascularization 
of the tumor and called "angiogenic capacity". The variation rate for 9 derived in [16] is : 

(3) 92{x,9) = cx - d9xi , 

If we denote X{t) = (x(t),9(t)) and define G(X) = (gi(x, 9), g2(x, 9)) we have the following system of 
ODE modeling the tumoral growth : 

(4) 

In the figure 1, we present some numerical simulations of the phase plan of the system. 
This system has been studied by A. d'Onofrio and A. Gandolfi in [10]. We define 




= G)' 



ft = (i,6) x (1,6), r = <9ft 



We will now turn our interest to the flow defined by the solutions of the system of ODE, as it will play 
a fundamental role in the sequel. We define the application 

$ . [o,oo[xr -> ft 

(t,<t) $ T (cr) 

as being the solution of the system (4) at time r with the initial condition a. We use of this application 
in order to see ft as ft ~ ]0, oo[xT. More precisely, we will show that <I> is an homeomorphism locally 
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Figure 1: Phase plan of the ODE system (4). 



bilipschitz ]0, cx)[xT* — > £1 where T* ;= T\{(b, b)}. In order to have a candidate for the inverse of <&, we 
define for (x,9) 6 

t(x,9) = inf{r > 0| $_ r (x,0) e V*} 
a(x,6) = <$>- T {xfi){x,9) 

The qualitative properties of the ODE imply the existence of such a couple (r, a) (the field points inward 
along r* and the solutions all converge to X* (see [10]) so going back in time they meet the boundary), 
and the Cauchy-Lipschitz theorem implies uniqueness because the system is autonomous and thus the 
characteristics don't cross each other in the phase plane. The time t(x, 9) is the time spent in O and 
a(x,9) is the entrance point of the characteristic passing through the point (x,9). From the Lipschitz 
regularity of we can't expect $ to be globally C 1 , this is why we introduce the following open sets : 

Oi = {$r(<x); ^r„r 6]0,oo[}, i = 1,2,3,4 

where 

r x =](i, i), (i, b% r 2 =](i, 6), (b, b)[, r 3 =](&, 6), (b, i)[, r 4 =](b, 1), (1, i)[ 

The restriction of $ to ]0, oofxTi is a diffeomorphism, as established in the following proposition. 
Proposition 2.1 (Properties of the flow). 

(i) The application $ is a diffeomorphism ]0, oo[xr.; — > Qi and for every r > and almost every a G T 

(5) Mr, a) = G ■ T^e-C div(G( *= W))ds 

where J$(r, a) is the Jacobian o/<f>. 

(ii) Globally, $ is an homeomorphism ]0,oo[xr* — > SI locally bilipschitz. 

Remark 1. The regularity proven here on $ validates the use of $ as a change of variables (see [12] for 
locally Lipschitz changes of variables). 

Proof. 

• <I> is one-to-one and onto. Let X = (x,9) £ f2. We have <5(t(Jl), c(X)) = X because &_ T /x)(X) = 
a(X) implies X = & T ( X ){°~{X)) (indeed <I>_ T is the inverse of 4> r when r is fixed). In the same way, 
(t($ t (<t)),o-($ t (o-))) = (t, ct). Thus $ is one-to-one and onto and <&~ l {x,9) = (t(x,9),o-(x,9)). 

• <fr is a diffeomorphism on ]0, oofxTj. Using the general theorem of dependency on the initial 
conditions for ODEs, $ is C 1 (]0, oofxTi) and if we call <r(s) a parametrization of Tj, we have fjf (t, er(s)) = 
Dj / ( I ) 7 -(ct(s)) oa'(s), and the following characterization of ct(s)) stands : for each s, it is the solution 
of the differential equation 

f f = DG($) o Z 
I ^(0) = 
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Using this characterization, we can derive the formula (5) for the Jacobian J$ (r, a) . We have J$ (r, a) 
d± 

dr 



¥ A f = ¥ A G ( $ )' and differentiating in r, we get 



9 <9$ a$ 

— Mt, a) = DG o — A G($) + ^ADGo G($) 
or as as 

= tracc(Z)G)J$(t,cr) = div (G)J^(t, a) 

Hence, for all <r(s), using that J*(0,ct(s)) = (t'(s) A G(a(s)) = |<r'(s)|G • T7(a(s)) ^ 0, we obtain the 
formula 



(6) 



J*(i,a(s)) = \<r'(s)\G-lt(a(s))exp( [ div(G($(r, a(s))))dr) ^0 

Jo 



We get (5) by choosing a paramctrization with velocity equal to one. 
Remark 2. In the sequel, we fix this parametrization 

We can then apply the global inversion theorem to conclude that $ is a C 1 -diffcomorphism ]0, oofxTj — > 
£V 

• Globally. From the given properties of the vector field G, we can extend the flow to a neighborhood 
V of f2, and we have that it is C 1 (]0, oo[xV) (see [8], XI p. 305). Hence $, which is the restriction of this 
application to ]0, oo[xr* with T* being Lipschitz, is locally Lipschitz. Remark here that it is not globally 
Lipschitz since ^$ T (cr) can blow up when r goes to infinity, due to the singularity at X*. 

To show that is also locally Lipschitz on £1 we consider some compact set K C £1 and show 
that <I> _1 is Lipschitz on K. We define Ki = £7; n K 7 and Ki :— <&~ 1 (K l ) c]0, oofxIV Now since $ 
is the restriction of a globally C 1 application, we have $ 6 C 1 (/f i ), meaning that its differential 1?$ is 
continuous until the boundary of Ki. Moreover using the formula (6), we see that the value of £>$ on dKi 
is invertible since we avoid the singularity X* . Hence, using the continuity of the inverse application we 
obtain that Z?$ _1 = (Z?^) -1 is continuous on Ki. Thus € C l (Ki) and so it is Lipschitz on each Ki. 
As the global continuity of on £1 is deduced from the continuity on £1 of X i-> t(X), it is Lipschitz 
on K. □ 



2.2 A renewal equation for the density of metastasis 

Starting from the velocity field G of the previous subsection for one given tumor, we now derive a 
renewal equation for the density p(t, x, 9) of metastasis at time t, size (= number of cells) x and so 
called "angiogenic capacity" 0. The term density for p means that the number of metastasis at time 
t in an infinitesimal volume centered in (x, 6) and of size dxdO is p(t, x,9)dxd0. We assume that each 
metastasis evolves in the space (x, 9) with the velocity G{x, 9). Expressing the conservation of the number 
of metastasis, we obtain 

(7) d tP + div(pG) = 0. 



The metastasis cannot have size nor angiogenic capacity bigger than the parameter b, and we assume 
them to have size and angiogenic capacity bigger than 1. We are thus driven to consider the transport 
equation (7) in the square £1 = (l,b) x (1,6). The field G pointing inward all along the boundary, we 
need now to precise the boundary condition on F. 

Metastasis do not only grow in size and angiogenic capacity, they are also able to emit new metastases. 
We denote by f3(x,9,a) the birth rate of new metastasis with size and angiogenic capacity a £ F by 
metastasis of size x and angiogenic capacity 9, and by f{t,a) the term corresponding to metastasis 
produced by the primary tumor. Expressing the equality between the entering flux of metastasis and the 
number of new born, we derive the following boundary condition on T : 

-G • l^(a)p(t, a) = I (3(x,6,o-)p(t,x,6)dxd6 + f(t,a) 
Jn 
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We then assume that there is no coupling between (x, 8) and a in the expression of (3, which is traduced 
by an expression of (3 as [3{x,6,a) = N(a)(3(x,9). Now let Q :=]0,+oo[xO, E :=]0,+oo[xr. The 
equation is 

( dtp + div(Gp) = V{t,x,6)€Q 

(8) I -G-l?p(t,<j) = N(a) f n /3{x,6)p(t,x,6)dxd6 + f(t,a) V(t,cr)e£ 

[ p(o,x,6) = p°{x,6) " v(M) e n 

We will do the following assumptions on the data : 

P £ L°°, P>0 a.e., N eLip(T) with compact support in T* , N > 0, J N = 1 

(9) G given by (2) and (3) 

Remark 3. In practice, the new metastasis only appear with size 1 and there should not exist metastasis 
on 1^2,3,4, thus in the biological model we have supp(iV) C T\. The expression of j3 we use in the 
biological applications is (3(x, 9) — mx a with m and a two positive parameters traducing respectively 
the aggressiveness of the cancer and the spatial organization of the vasculature. The source term f has 
the following expression in biological applications : f(t,a) = N(a)f3(X p (t)) with X p (t) representing the 
primary tumor and being solution of the system (4) . 

Definition 2.1 (Weak solution). Let p° e L 1 (il) and f e L 1 (]0, oo[xT). We call weak solution of 
the equation (8) any function p G C([0, co[; L 1 (il)) which verifies: for every T > and every function 
4> e C^([0,+oo[xn*) 

(10) f [ p[d t <f> + G-V<f>]dtdxd9+ [ p°{x,9)(t>{Q,x,9)dxd9 
Jo Jn Jn 

p{T 7 x, 9)<j)(T, x, 9)dxd9 - [ [ N(a) ( [ f3(x, 6)p(t, x, ff)dxdJd \ <f>(t, o)dodt = 
! Jo Jr \Jn J 

Analyzing the equation (8) indicates that the solution is the sum of two terms : an homogeneous one 
associated to the initial condition, which solves the equation without the source term / (which we will 
refer to as the homogeneous equation) 

f dtp + div(Gp) = V(t,M) G Q 

(11) <^ -G • l?p(t, a) = N(a) J n (3p(t)dxd9 V(t,cr)eS 
[ p(0,x,9) = p°{x,9) ' "" V(x,0) G Q 

and a non-homogeneous term associated to the contribution of the source term f(t,a) and solution to 
the equation (which will be refered as the non-homogeneous equation) 

f dtp + div(Gp) = V(t,M) e Q 

(12) <^ -G-l^p(t,a)=N(o-)J n pp(t)dxde + f(t,a) V(t,cr)eE 
[ p(0,x,0) = V(x,9) en 

For existence and uniqueness of solutions, we will deal with the homogeneous problem using the semigroup 
theory and with the non-homogeneous one via a fixed point argument. 

2.3 Semigroup formulation for the homogeneous problem 

We reformulate (11) as a Cauchy problem 



(13) 



d tP {t) = Ap(t) 
p(0) = p° 
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We introduce the following space: 



Wdiv(ft) = {V E div(Gy) £ L 1 ^)}, 

and the following operator 

V i ^ -div(GV) ' 

where 

(14) D(A) = {V G W div (n); -G.lt ■ j(V)(a) = N(a) [ (3{x, 9)V(x, 9)dxd9 1 Vcr G T} 

JO 

We refer to the appendix for a short study of the space Wdiv(^) 5 in particular the definition of the 
application ~/(V) as the trace application. 

There are three definitions of solutions : the classical (or regular) solutions, the mild solutions ([15] 
II. 6, p. 145) and the distributional solutions (2.1 with the source term / = 0), the second and third ones 
being two a priori different types of weak solutions. The following proposition proves that the weak and 
mild solutions are the same ones. 

Proposition 2.2. Let p G C([0, oo[; L 1 (Q)), then 

(p is a mild solution of (11)) (p is a weak solution of (11)) 

Proof. • First implication =>- : It comes from the fact that mild solutions are limit of classical ones which 
are weak solutions in the sense of definition 2.1, by passing to the limit in the identity (10). 
• Second implication <= : Let p G C([0, oo[; L 1 (0)) be a weak solution in the sense of definition 2.1 with 
/ = 0. Define the function R{t) = J Q p(s)ds. We verify now that R(t) G Wdi V (^) by using the definition. 
Fix t > and a function <f> E C^(Q). Using the function <fi(t, x, 9) = (f>(x, 9) in (10), we have 

p(s)ds(G ■ V<t>)dxd6 = - [ p°{x,9)(f>{x,9)dxd9 + [ p(t,x,9)4>{x,9)dxd9 

Jn Jn 

Therefore R(t) G W div (0) and p(t) = A f* p(s)ds + p°. 

We now prove the boundary condition part contained in order to have R(t) G D(A). Let <p(a) be 
a continuous function on F, with compact support in T*. We can extend it to a function of C c (f2 ), 
still denoted by <f>, by following the characteristics and truncating, namely : <p(^ T (a)) — 0(ct)C(t), t G 
[0, +oo[, a G T with ((t) being any regular function with compact support in [0, +oo[ such that £(0) = 1. 

Now, using the density of C^(0*) in C C (Q*), choose a family <p e G C^(0*) such that (j> e — — > <f>. For each 
e, using the remark following the definition of weak solutions with the test function (f> e (t, x, 9) = 4> e (x, 9), 
we have for every t > 

f R(t)G-V(f) t + [ p°{x,6)4> € {x,6)dxd6- [ p(t,x,9)cf> t (t,x,9)dxd9 = 
Jn Jn Jn 

I N{a)(j) e {a)da [ f3(x,9)R{t)dxd9 
Jr Jn 

As R(t) G Wdiv(^), and — div(Gi?) = p — p° by passing to the limit in e, we obtain 

/ ~/{R(t))G -lt(f>= I N<p I /3R, Vt > 
Jr Jr Jn 

This identity being true for any function <j> G C C (T*), we have the required boundary condition on R(t). 
This ends the proof. □ 
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3 Properties of the operator 



We first remark that (A,D(A)) is closed, by classical considerations and the continuity of the trace 
application (prop.A.l). 

3.1 Density of D(A) in L\Q) and adjoint (A*, D(A*)) of the operator 
Proposition 3.1. The space D{A) is dense in L x {£i) 

Proof. The proof follows the one done in [2] in dimension 1, although some technical difficulties appear 
in dimension 2. Since C x (fl) is dense in it is sufficient to approximate any function / G C x (fl) by 

functions of D(A), for the L 1 norm. Thus let / E Cl(tt) be a fixed function. Let E CC T* = T\(b, b) be 
the support of N(o-) and for each n E N let V n be an open neighborhood of £ such that mes(V n ) — > 0, 
and (6, 6) ^ V n . There exists a function <fr n E C^(M 2 ) such that 

*.<*.*>={; if ( L%)^ 

Then, we extend the function H(a) = _q.^^ ■ r* ^ R to a Lipschitz function H : V n fl fl — > K (for 
example by following the characteristics). Let 



h n {x,6) 



{H<f» n )(x,e) tf(x,6)e\^nn 
o if (x, e)n\v n 



It satisfies /i„ € W 1,00 (Q) and /i„ — » 0. Let /„ = / + (i„/i n , with 



1 - J n j3h n dxd0' 

Since H/inHi^m) — ^ and /3 is in L°°, for n sufficiently large 1 — | J n /3h n dxd0\ > 1/2 and \a n \ < 
2||/3||loo||/||li. Then /„ ^ f and furthermore, since h n E W 1 ' 00 ^) C W 1 * 1 ^) C W div (0), we have 

/n G D(A). □ 

We are now interested in characterizing the adjoint of the operator (A, D{A)). We will see that the first 
eigenvector of (A*,D(A*)) plays an important role in the structure of the equation in the asymptotic 
behavior (see theorem 4.1). 

Proposition 3.2 (Domain and expression of A*). 

(15) D(A*) = {UeL°°; G-VU E L°°} := Wj? (fi) 

A*U = G-X7U + I3 J U(a)N(a)dcr. 

Proof. • The first inclusion for the domain of A* is a consequence of the property A.l. The second 
inclusion D(A*) C W^ v (Cl) requires a little much of work. For a function U E D(A*), we will show that 
4> M>< U, div(G<f)) > can be extended in a continuous linear form on L 1 , which will allow us to conclude 
using the Riesz theorem that U E W£? v . To do this, it is sufficient to show that there exists a constant 
c > such that 

(16) \<U,A(j) > V ',v | < c\\<j)\\ L i V<f>£V(fl) 

This is almost done by the definition of the domain D(A*) except the fact that T>(Q) is not a subset of 
D{A). We are driven to use the following trick. Define the space : 



T> (fl) = {</>€ T>(fl); / /30 = 0} 
Jn 
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which is a subspace of D(A). We will project a given function in D(0) on T>o(Cl). Let <\>\ e P(0) be a 
fixed function such that / fi(f>i = 1. Then 



n - 



£V (n)cD(A) £R<fii 

So eventually, denoting as c\ the constant given by the belonging of U to D{A*~) 

< U,A<\> > v , tV | = | < U,A(<b-(^J Pc^j 0i) > + < ?7, (y*/?^ A0i > | 
<(ci +c 1 ||/3|| i oo||0 1 || il + ||/3|Uoo||C/|| i oo||^ 1 || il )|H| L1 
which shows (16) and thus yields the result. □ 

3.2 Spectral properties and dissipativity 

In order to have a candidate for a stable asymptotic distribution of the solutions of our equation, we are 
interested in the stationary eigenvalue problem : 

f (A, V, $) e I| x D(A) x D(A*) 

(17) \ AV = XV, A*tf = Atf 

[ J n V^dxdO = 1, * > 0, / r NVdtr = 1 

Proposition 3.3. [Existence of solutions to the eigenproblem] Under the assumption 

(18) f f l3(^ T {a))N(a)dTda > 1, 

Jo Jr 

there exists a unique solution (Ao,V, *) to the eigenproblem (17). Moreover, we have the following 
spectral equation on A : 



(19) 



i> + oo p 

J J P(^ T (a))N(a)e- XaT dTda = 1 



The direct eienvector is given by 

V(<$> T (a)) = Cx N(<7)e- XoT \J^\-\ Vt > 0, a.e a E V 

where C\ is a positive constant and |J$| is the jacobian of $ from section 2.1. The adjoint eigenvector 
^ is given by : 

/oo 
/3(^ s (a))e- XoS ds Vt > 0, a.e a e V. 

Hence we have 

-r , r\ \ sup/3 , . ^ 
— — < *(z, 6>) < — — V(x, 9) eil 
Ao Ao 

Remark 4. In the model we use in practice, where f3(x, 9) = mx a the condition (18) is fulfilled since 

/ / f3(^ T (a))N(a)dTd(T = oo, and the inequalities on ^ write 
Jo Jr 

t- <*(x,6) <— V(x 1 6)en 
Ao Aq 
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Proof. • The direct eigenproblem. 

o We use the following change of variable, which consists in transforming a function of Wdi v (^) into a 
function of W 1 ' 1 ^, +00) ; L 1 (r)) : 

V(t,*) = -V(* t (*))\J 9 \, Vtg[0,+oo[, aeT 

where we recall that |J$| = — G ■ T^(cr)e^o dlv (G(<M<r)))ds j s ^he j aco bian of the application $ : (r, a) 
3>t(c) (see section 2.1). 

Rewriting the problem on V and denoting /3(r, cr) = /3($ T (cr)), we get 



<9 T U + AU = 
1/(0, cr) = Ar(cr) J pVdrda' 



(21) 

Direct computations show that Problem 21 has a solution if 

(22) l = jf J N ( (T )P( T > cr ) e ~ XTd < jdT 

and conversely, if A is a solution of the equation (22), we get solutions to the problem (21) given by 

(23) V(r,a) = C Xo N(a)e- x ° T 

and we can then fix the constant C\ a > in order to have the normalization condition 1 = J n V^dxdQ 
with \& the dual eigenvector defined below. 

o We now prove that there exists a unique solution to equation (22) under the hypothesis (18). Indeed, 
let us define the function F : M -> R by 

(24) F(\)=J (^J N(a)P(T,afj e'^dadr 

It is the Laplace transform of the function t J t N(a)j3(r, a)da. The condition (18) means that 
F(0) > 1 and F being strictly decreasing on K and continue on ]0, +oo[, the equation (22) has a unique 
solution in M, A s]0, +00 [. 

o From (23), we obtain that V <E T4 /1,1 ((0, +00) ; L 1 (r)). Using theorem A.l from the appendix, we 
deduce that V £ W d w(ty- 

Remark 5. Here the theorem A.l takes its interest since it is not completely obvious that the composition 
of V by would give a function in W / div(^), due to the fact that the change of variable $ (and 
is not globally Lipschitz. 

• The adjoint eigenproblem. 

o Expression of 'P. Using the expression of the adjoint operator (A* , D(A*)) from the proposition 3.2, 
the adjoint spectral problem reads, along the characteristics : find G (Q) such that 

(25) a T *($ T (cr)) = A *($ T (cr)) - /3($ r (a)) J * {a')N{a')da' 

from which we get, for each function ^f(a) defined on the boundary, a solution to the equation given by 

(26) f($ T (<7)) = *(cr)e AoT - f ^{a')N(a')da' f P(^ s (a))eM T - s '> 

Jr Jo 

o Non-negative solution. To get a non-negative solution we are driven to the following condition 

tf(tr) > f ^{a')N{a')da' [ f3(^ s (<r))e- x " s ds, a.ecreT 
Jr Ja 
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Now, if the inequality is strict, multiplying by N(a) and integrating on T gives 1 > / r J °° f3($ s (<j))e x ° s dsd<r 
which belies the spectral equation (22). We are thus driven to choose 

(27) *(<7)= / ^{a')N(a')da' [ P(<S> s (o-))e~ x ° s ds, Ver 6 T 

Jr Jo 

Defining g(a) = J* °° /3(^ s (a))e- x ° s ds, this means that ^(cr) is in the vector space generated by g. Then 
it remains to have the suitable normalization constant. Remembering the spectral equation (22) verified 
by Ao shows that the function = g(a) is appropriate. We finally get (20) from (26), which gives 

*eL°° and ||*||l"o < ^t^- 

o Regularity of Using the equation (25) verified by \& we get d T ^f(^ T (a)) G L°° and so using the 
conjugation theorem of W£> V (Q) and TU 1 ' oo ((0, oo) ; L°°(T)) (theorem A.l), we have * e W|? (Q). □ 

Using the change of variables V(t, <j) — V(& T (<r))| J$|, the theorem A.l and the proposition A.l, we can 
follow the methods of the one-dimensional case done in [3, 2] thanks to the decoupling of (3(x,6,a) in 
N(a) x (3(x,0), to obtain the following proposition. 

Proposition 3.4. (i) For Re(X) > Ao, we have Im(XI — A) = L 1 ^!). 
(ii) The operator (A — uI,D(A)) is dissipative for every uj > 

Applying the Lumer-Philips theorem, we obtain 

Corollary 3.1. Under the assumptions (9) the operator (A,D(A)) generates a semigroup on L 1 (f2) 
denoted by e tA and we have 

\\\e tA \\\ < e^^ 00 



4 Existence and asymptotic behavior 

4.1 Well-posedness of the equation 

4.1.1 Existence for the non- homogeneous problem 

Proposition 4.1. 

(i) Let f E L 1 (]0, oo[; i x (r)) and assume (9). There exists a unique solution of the non-homogeneous 
problem (12), denoted by Tf and we have 

Tf eC([0,oo[;L\Q)) 

(ii) If f € ^([0,001; L 1 ^)) and /(0) = then 

Tf e C\[0, oo[; L 1 ^)) n C([0, oo[; W div (0)) 
Moreover, we have the positivity property 

(/ > 0) =► (Tf > 0) 

Proof. The proof is based on a fixed point argument. It is divided in three steps : first we prove the point 
(ii) using the Banach fixed point theorem, then thanks to an estimate in C([0, oo[, L 1 (ri)) we construct 
the weak solutions as limits of regular solutions, and finally we prove uniqueness. 

• Step 1. o As usual now, we first simplify the problem using the conjugation theorem (theorem A.l). 
We use the change of variable p — p(Q T (o-))\ J$ and still denoting p for p and f3 for /3 = /3($ T (er)), we 
consider the following non-homogeneous problem with nonzero initial condition 

{d t p + d T p = 
p(t,a) = N(a)Jpw + f(t,a) 
p(0) = p° 
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Let p° G D(A) and / G £^([0, oo[; L 1 ^)) with /(0) = 0. For T > we define the space 

X T = {w G C 1 ([0,T];L 1 (]0,oo[;L 1 (r)); w(0, •) = p } 
It is a complete metric space. To w G Xr we associate the solution p of the equation (28), namely 

p(t, t, a) = |iV( CT ) - r, r', crVr'rfa' + f(t - r, <r) |l t>T + P °(t - t, a)l t<T 
. ' 

—H(t-T,a) 

and define the linear operator T p oj by T p »jw := p. Note here that w > implies p > if p° > and 
/> 0, and that H G ^([O, T]; ^(T)). 

o Regularity of p. We now show that p E X T and that p G C([0, T]; IT M ((0, +oo) ; L 1 ^))). Indeed we 
have 

(29) p(t,T,a)l t>T = H(t-T,a)l t>T , p(t, r, a)l t<T = p°(r - t, a)l t<T . 

From these expressions we get that p G C([0, T]; L 1 (]0, co[; L 1 (r))) since the two functions H and p° are 
inL 1 . 

Moreover, H(0, a) = N(a) J (3p° = p°(0) from the compatibility conditions contained in the facts that 
w G X T , /(0) = and p° G D(A). This allows to conclude that p(t,-) G C([0, oo[; L 1 ^)). Further- 
more, from the expressions (29), we see that for each i p(i, •) G W 1 ' 1 ^, t), L 1 ^)) n IT 14 ((t, oo), L x (r)) 
since p° G -D(^l) and H G C 1 ([0, T]; i 1 (r)). Combining this with the continuity in r gives p(t, •) G 
VF 1,1 ((0, +oo) ; L x (r)). Finally from the expression of d T p obtained differentiating in r the expressions 
(29) we get p G C([0, T], ^((O, +oo) ; L 1 ^))). 

It remains to show that p G C 1 ([0, T]; L 1 (]0, oo[; L x (r))). For the sake of simplicity we forget the 
dependency on a. We define for almost every t and r 

ftp(t,T) := d t H(t - T )l t>T - d t p°(T - i)l t<T 

+ - /l9tp(t)||l,i([0,oo[) = 

±||ff(t + h - •) - tf(t - •) - hd t H(t -)IUh[o,*[) 

+ + h - •) - p°(- - t) - ^p°(r - i)|| Ll(]t , t+M) 

S v ' 

A 

+ l\\P°(- -t-h)- p°(- - t) - ^p°(- - t)\\ L i {[t+h>oo[) 

The first and the last terms go to zero when h tends to zero since H is in C 1 ([0, T]; L 1 (]0, oo[)) and p° is 
in D(A). To deal with the last term A, we write 

A ^J l \ t \H(t + h-T)- P °(T-t)\dT + J d tP °(T-t)dT 

The first term goes to zero because of the compatibility condition H(0) = p°(0) and also the last one 
because d t p° G L 1 . We can then conclude p G C 1 ([0, T]; L 1 (]0, oo[)). 

o The previous considerations show that the operator 7^,0 / has values in Xt- Now, if wi and w 2 are in 
X T we compute 

\\TpOjW! - Tp0jW 2 \\x T < T \\/3\\l^\\wi - W 2 \\x T 

Using a bootstrap argument we prove the existence of a solution on [0, oo[ and transporting the regularity 
facts back to ft by using the conjugation theorem A.l ends the point (ii). 

• Step 2. Denote by Tf the fixed point of the operator Toj, defined up to now only when / is regular 
and satisfies the compatibility condition /(0) = 0, one has 



Now we compute 
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Lemma 4.1. Let f E C 1 ([0, oo[; L 1 (r)) with /(O) = and Tf be the solution of the equation (28) with 
a zero initial condition. Then for all T > 



l|T/|| c([ o.T]:Li(n)) <e T "^ f \f{s)\e-WW~ s ds 

Jo 



Proof. The solution Tf = p being regular, the function \p\ also verifies the equation and integrating on 
f2 yields 



d_ 

It 



[ \p\(t)drda=\ [ !3p(t)dxd9+ [ f(t, a)da\ < f \p\(t)dxd9 + f \f(t,a)\da 

Jn Jn Jr Jn Jr 



and a Gronwall lemma gives the result. □ 

Now by a density argument and using the previous lemma, we can construct a solution Tf E C([0, oo[; L 1 (ft)) 
when / G L 1 (]0 , oo[; L 1 (T)) . This solution can be constructed non-negative whenever / is non-negative 
itself. 

• Step 3. It remains to show the uniqueness of the solution. If p\ and p 2 are two solutions of the non- 
homogeneous equation (28), then p\ — p 2 is a weak solution of the homogeneous equation (11) with zero 
initial condition. From the proposition 2.2 the weak solutions in the sense of the distributions are the 
same than the mild solutions and thus pi — p 2 is a mild solution of the homogeneous equation and hence 
is zero by uniqueness of the mild solutions □ 

4.1.2 Existence for the global problem 
Theorem 4.1 (Existence and uniqueness). 

• Let p° E L X (Q) and f E L 1 (]0,oo[xT), and assume (9). There exists a unique weak solution of the 
equation (8), given by 

p = e tA p° + Tf 

with Tf being a weak solution of the non-homogeneous equation (12). 

• If p° E D(A) and f E C 1 ([0, oo[; i 1 ^)) and verifies /(0) = 0, then we have 

p e c\[o, oo[ ; LHn)) n c([o, oo[ ; w div (n)) 
4.2 Properties of the solutions and asymptotic behavior 

In the next proposition, we prove some useful properties of the solutions, which appear in the L\, norm 
defined by 

(30) 11/11^=^1/1*^, 

with * the dual eigenvector from proposition 3.3. We should notice that when (3 E L°° and (3 > 5 > 0, 
by the inequalities from proposition 3.3, the L\, norm is equivalent to the L 1 norm. Hence the solutions 
have finite L\, norm. The main idea in the proof of the following proposition is to use various entropies 
in the space L\, and is based on ideas from [20] and [19] 

Proposition 4.2. Let p a E i 1 (fi) and p the solution of the equation (8). The following properties hold : 
(i) 

(31) / \p(t)\H> <e A °'( / /* / *(o-)e- x ° s \f\(s,o-)do-ds\ , Vt > 

Jn Un Jo Jr J 
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(ii) (Evolution of the mean-value in L\,) 

f p(t)^ = e Xot \( p°V + f f ^{a)e- XoS f( Sl a)dads\ , Vi > 
Jo Un Jo Jr J 

(Hi) (Comparison principle) If f > 

P?<P2 Pl(t)<P2(t) Vt>0 

Proof. Each time we aim to prove something on weak solutions, we will start proving it for classical 
solutions and then use the density of D(A) to conclude. So, let do the calculations with a strong solution 
p associated to an initial condition p° in D(A) and a function / G C 1 (]0, oo[; L 1 (r)) with /(0) = 0, for 
which the calculations can be justified. We first remark that the dual eigenvector ^ which belongs to 
W£? V (Q) verifies the following equation : 

(32) G • V* - A * = -,3 

since by the construction of \f and the spectral equation J r ^(cr)N(a)da = 1. Defining p(t,x,6) = 
e~ Xot p(t, x, 9) we have the following equation on p : 

(33) d t p + div(Gp) + \ p = 0, 

with the same initial condition as for p and a suitable boundary condition. Using that p G Wdiv(^l), 
* G W^ v (Sl) and the proposition A. 2, we obtain the following equation on p^ : 

(34) d t {p^) + div(Gp*) = -0p 

(i) Let us first state the following lemma. 

Lemma 4.2. Let p° G L 1 ^) and p be the associated weak solution of the equation (8) . Then the function 
\p\ solves the same equation, with suitable initial and boundary conditions. 

Proof. For a regular solution of the equation p associated to a regular initial condition p° G D{A) and a 
regular data /, we can use the proposition A. 2 with the function H(-) = | • | to have that \p(t)\ G Wdi V (^) 
and 

div(G|p|) = sgn(p)G • Vp + |p|div(G) 

Since p is regular in time, by multiplying the equation by sgn(p) we get the result. For a solution 
p{t) G L 1 (r2) we obtain the result by density of the strong solutions. □ 

Thanks to this lemma we have the equation (34) written on \p\, from which we get, integrating in (x, 6), 
that 

^ f \p\»fdxd6 = - [ j(\p\)^G-~^dcr- [ P(x,6)\p(t,x,6)\dxdB 
dt Jq J r Jn 

= [ *(tr) N(a) [ f3(x,6)p(t,x,6)dxd9 + e- Xot f(t,cr) - [ 0(x,6)\p\(t,x,6)\dxd0 
Jr Jn Jn 

<£\f(t,a)M*) 

and thus deduce the first property by integrating in time. To deal with weak solutions we again use the 
density of regular solutions. 

(ii) To obtain the evolution of the mean value, we integrate in space and use again a density argument. 
(Hi) Writing the solution of the global problem as p — e tA p° + Tf, we only have to prove the positivity 
for the homogeneous part since the positivity of the non-homogeneous one has been established in the 
proposition 4.1. It can be proved in the same way as the first point but using the negative part function 
instead of the absolute value. □ 
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Proposition 4.3 (Asymptotic behavior). Assume that 

J /3($ T (a))N(a)dTda > 1, 

and that there exists p > such that (3 — p^ > 0. Let p° e L 1 ^), f e L 1 (]0, oo[xr), p the associated 
solution to the global problem and (Ao, V, ^) £ x D(A) x D(A*) be solutions to the direct and adjoint 
eigenproblems. We have : 



\\p(t)e 



-A t 



- m(t)V\\ L i < e""* - moV\\ L i + 2 J* e~^-^' jf |/|(«, a)tt(<r)<k 



where \\f\\ L i = f and m(t) = f p(t)V = f p ^ + f e- x ° s [ f(s,a)^{a)dads. 

Remark 6. Notice that choosing p < Ao swes £/ie convergence of the integral J °° e~( Ao_ ^ s J r |/|(s, <7)\P(ct)</s 
and i/ius i/ie convergence to zero of the right hand side of the inequality. 

Remark 7. The hypothesis of the theorem are fulfilled in the case of biological applications where 
P(x,6) = mx a , because we have then f3 > m > and ^ £ L°° . 

Proof. Again we start with a regular solution p(t, x, 9). We then follow the calculation done in [20] III. 7 
pp. 66-67, adapting the method to take into account the contribution of the source term. Define the 
function 

h(t, x, 9) = p(t, x, 9)e- Xot - m(t)V 

which satisfies J a h(t)^ = for all non-negative t, by the property of evolution of the mean value and 
since J Q = 1. As the direct eigenvector V solves the equation (33), h solves the equation 

d t h + div(hG) + \ Q h = - e - Xot FV 

where F(t) := J f(t,a)^>(a). Multiplying the equation by the function sgn(/i) gives the following 
equation on \h\ 

t \h\ + div(|/i|G) + Ao|/i| = -e- Xot FVsgn(h) 
Multiplying this equation by ^P, the equation on <]/ by \h\ and then summing the both gives 

dt(\h\V) + div(G|fe|tf) = -f3\h\ - e - Xot FV^sgn(h) 

Now integrating in (x, 9) yields : 



Now we use that 



^ / \h\*dxd0 = [ N(a) [ (3hdxd9 + e- Xat f(t,a) 

dt Jn Jr Jn 

-e- Xat F I V^sgn(h)dxd9 
Jn 



da 



(3\h\dxd9 



N(a) I (3hdxd9 + e- Xot f(t,cr) 



= (n(ct) J Phdxd9^j sgn(h(cr))+ 
(e-^ffra)) S gn(h(a)) 



to obtain 



d_ 

~dl. 



[ \h\^dxd9 = [ «f(a)N(a)sgn(h(a))do- [ (3hdxd9 
Jn Jr Jn 



(3\h\dxd9 



+ [ e- Xot f{t 7 (j)sgn(h(o-))^(a)do--e- Xot F [ V^sgn{h)dxd9 
Jr Jn 
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We first deal with the term A. Using that J Q = and remembering that J T Nty = 1 we compute 



A < 



[ Phdxd6-n [ ^hdxdO - [ P\h\dxd9 < [ {(3 - y$>)\h\dxd6 - [ (3\h\dxd9 
Jn Jn Jn Jn Jn 



< -At / ^\h\dxd6 
Jn 

where we used that (3 — \i§> > 0. 

• A direct majoration, the positivity of the eigenvectors V and \f/ and the fact that / V 1 ® = 1 gives, 

Jn 



denoting F(t) := J \f(t, cr)|*(a), that B < 2e- Xot F(t) 
• A Gronwall lemma finally gives 



/ \h(t)\V < e-"* { [ |/i(0)|* + 2 / e- (x °-^ s F(s)ds 
Jn [Jn Jo 



which is the required result. For an initial data in remark that it is possible to pass to the limit 

in the previous expression. □ 

5 Conclusion and perspectives 

In the present paper, we achieved the first step of our program consisting in elaborating and applying 
a model of metastatic growth including the tumoral angiogencsis process : the mathematical analysis 
of the direct problem. To do this, we used semigroup techniques and also the characteristics in order 
to study the natural regularity of the solutions to our equation, which led us to a short study of the 
space Wdiv(^)- This theoretical study brings to light the quantity Ao as characterizing the asymptotic 
growth of the metastatic process. This parameter has biological relevance and finding the best way of 
controlling its value by means of antiangiogenic drugs can be of great interest. The crucial problem is now 
the identification of the parameters of the model from biological data, in order to predict the optimized 
administration protocol for antiangiogenic drugs. 



To achieve this, we need to perform efficient numerical simulations of the equation. Due to the large 
disproportion of the boundary condition and the solution itself, as well as the size of the domain (typically 
6 = 10 11 for humans) and the behavior of the characteristics attached to the velocity field G (see figure 
1), performing good simulations of the equation is not an easy task. In particular, classical upwind 
schemes are not efficient. We are currently working on a characteristic scheme which follows the one 
used in [3]. We will then include the anti-angiogenic treatment in the equation, which mathematically 
means transforming G in a non-autonomous vector field. We will also address the inverse problem and 
the parameter identification. Our model has the good property that it has a small number of parameters. 
So we hope that it can be used efficiently to make predictions. We want to study mathematically the 
parameter identification. 

A A short study of iy d i v (^) 

Let 

n = (1,6) x (1,6), r = on 

and G the vector field on ft with components given by (2) and (3). Consider the space 

W div (n) := {V G L\n) | 3g e L^fi) s. t. f VG ■ V<£ = - / g<f>, V<f> G C^(fi)} 

Jn Jn 
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The function g of this definition is denoted div(GT^). We endow this space with the norm 

\\V\\ w ^ = \\V\\ L , + \\div{GV)\\ L , 

With this norm, Wdi V (^) is a Banach space. In the following we also denote this space by W^ iv (f2). 

Remark 8. //div(G) G L°° and V G Wdi V (^), we can define 

G ■ W := div(GV) - Vdiv(G) G L\n) 

and the space Wdw(Q) is also the space of L functions such that there exists a function g G L (fl) 
verifying 

f vdiv(G4>) = - [ g <t>, y^ecKn) 

Jn Jn 

This space already appeared for the study of the boundary problem for the transport equation (see 
[4, 6, 7]). 

In the same way, we define the space 

w^ v (n) = {u g l°°(Q)| g-vu e L°°(n)} 

A.l Conjugation of W div (Q) and W 1 ' 1 ^, +oo) ; L 1 ^)) 

For a function V G L 1 , the fact of belonging to Wdiv(^) means that it is weakly derivable along the 
characteristics. The next theorem makes this more precise. 

Theorem A.l (Conjugation of Wdi V (^) and W 1 ' 1 ((0, +oo) ; L 1 (r))). Let p = 1 or oo. The spaces 
Wj iv (0) and T4 /1,p ((0, +oo); L P (T)) are conjugated via $ in the following sense : 

V G W^ w (fl) & (Vo*)|J»|Vp e ^((o.+oo); LP(T)) 

Moreover, for V G W^ iv (fl) we have almost everywhere 

d T (Vo<p\j^/p) 

and the application 



(div(GV)o$)|J$| ifp=l 
(G • W) o $ if p = oo 



W£ v (fi) -> ^((0, +oo); L"(T)) 
V i y Vo*|J»|Vp 



is an isometry. 



Remark 9. 

• In particular, we deduce from the theorem applied to the function V = 1 i/iat | J$| G M /1 ' 1 ((0, +oo) ; i 1 (r)) 
and we recognize the well known formula 

dr\J$\ = div(G)|J$|, a.e 

• Smce by proposition 2.1, we have |J$| _1 G W^^QO, oo[xr*) we deduce that 

V(* T (v)) = n$ T ((7))|J*| x iJ^r 1 G Wi^ao,oo[; LL(r*)) 

d T V(* T (v)) =G-VC/($ T (a)) 
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Proof. We first show the theorem on W div (n) and then for (Q) 

• We prove now (V e W div (0)) (V := (V o $)|J$| e ^((O.+oo) ; L 1 ^))). Let V £ W div (Q) 
and remark that V £ ^ 1 (]0, oo[xr) since J$| is the Jacobian of the change of variable between f2 and 
]0, oo[xT*. Then, using the definition of VF 1,:L ((0, +00) ; i 1 (r)) we have to prove that there exists a 
function g £ L 1 (]0, oo[xT) such that for every function -0 £ C£°(]0, 00 [) 



/ V(r,a)^'(T)dT = - g(r,a)^( 
Jo Jo 



r)dr, a.e. a £ T. 

As we aim to use the change of variable $ which lives in ]0, 00 [xT, we will rather prove that for every 
function ( £ LipiT) (the Lipschitz functions on F) 

(35) jf |^°° V(t, a)^(r)drj ((a)da = 

J {-J div(GV)(* T (cr))|J$|^(r)dr| ((a)dcr 

which is sufficient to prove the result. Let now define the function 

4>(x,9) :=$(t{x,0)) 

with t(x,6) the time spent in ft defined in the section 2.1. Then ?/> has compact support in and 
is Lipschitz as the composition of a regular function and a Lipschitz function (see prop. 2.1 for the 
locally Lipschitz regularity of the function (x,0) M> t(x,9)), thus diffcrentiable almost everywhere and 
the reverse formula ip(r) = 4>(<fr T (l, 1)) (or ip(<5> T (a)) for any a £ V since the function ip depends only on 
the time spent in Q) yields 

V>'(t) = G ■ V^($ T (1, 1)), a.e. re]0,oo[ 
since ri->$ T (l,l)isC 1 . Doing now the change of variables in the left hand side of (35) yields 

(36) f\f V(T,a)^'(T)dT\c(o-)da = f V{x,6)C(o-(x,6))G ■ Vip(x,6)dxd6 
Jr I Jo J in 

Still denoting £(x, 6*) the function C( cr ( a; ! we remark that this function only depends on the entrance 
point a(x, 9) and thus we have 

(G • VC)($») = 9r(C(*r W)) = 9r(CW) =0, Vr > 0, a.e a 

To pursue the calculation, we need to regularize the Lipschitz functions £ and ^ in order to use them in 
the distributional definition of div(GV). We use the following lemma, whose proof can be found in [25], 
p.60. 

Lemma A.l. Let f £ W 1,0 °(fl) with a Lipschitz domain. Then there exists a sequence f n £ C°°(fi) 
such that 

w 1,p 

fn > f VI < p < oo, /„->•/ L°° weak - *, V/„ -> V/ L°° wea/c - * 

Now let ip n — > ip and Cm C as in the lemma. From the demonstration of the lemma which is done by 
convolution with a mollificr, since ip has compact support, so does ip n for n large enough. Now remark 
that for each n and m 

G • V(ip n ( m ) = CmG • W>„ + V«G • VCn 
The function VnCm i s now valid in the distributional definition of div(GV) and we have 

f V( m G ■ V^p n dxd9 = [ VG ■ V(^ n ( m )dxd9 - [ Vip n G ■ V( m 
Jn Jn Jn 

= - [ div {GV)ip n ( m dxd6 - [ Vip n G ■ \7( m dxd9 
Jn Jn 
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Letting first n going to infinity, then m and remembering that G • V( = yields 

/ V(G ■ VtpdxdO = - [ 6xv{GV)i)C,dxd6 
Jn Jn 

Now doing back the change of variables gives the identity (35). 

• We show now the reverse implication. Let V £ VF 1 ^ 1 ((0, +00) ; i 1 (r)) and tp £ C£°(Q). Define 
V(x, 6) := (V o $ _1 )|J $ -i| and tp(r,a) :— tp(^ T (a)). Hence tp is C\ in the variable r and we have 
d T if = (G • W) o $. Now 

VG-Vipdxd0 = / V(T,a)d T xp{T,a)dTda = / d T Vipdrda 

Jn Jr J a Jr Jo 

= - [ d T V o $- 1 il>dxd6 
Jn 

Hence we have proved that V £ Wdiv(^) and that div(GV) = d T V o 

• We prove now the part of the theorem on WJ? (fi). Let U £ Wj? (fi) C W div (Q). Then £/ o $ £ 
L°°(]0, 00 [xL). Moreover, following the second point of the remark following the theorem, we have 

d T U(<f> T (a)) = G ■ VE/($ r (cr)) £ L°°(]0, oo[xL). 

Using that for U £ W 1>oo ((0, +00); L°°(T)) we have locally G-VU := o^t/o^- 1 e with [/ = C/o*- 1 

gives the reverse implication. □ 

A. 2 Trace theorem, integration by part and calculus of functions in Wdi V (0>) 

Thanks to the theorem A.l, we can now transport the theory of vector- valued Sobolev spaces to W-div(^), 
for which we refer to [12]. 

Proposition A.l (Trace in Wdiv(^) and integration by part). Let p = 1 or 00 and V £ W^ iv (Cl). We 
call trace of V the following function 

■y(V)(a) = (Vo<f>)(0,a), Vcr £ T 

We have / y(V)G ■ ~T? £ L P (T) and there exists C > such that 

\h(v)G-^\\ LP(r) < c\\v\\ wSiv{n) , w £ wsjn) 

Moreover, if V £ W diw {Q) and U £ W^ V (Q). Then 

[ [ Udiv(GV) + [ [ VG-VU = - f 7(V)7(I7)G • ~& 
J Jn J Jn Jr 

Proof. It is a direct consequence of the properties of functions in W 14 ((0,+oo) ; L 1 ^)) and in ^'^((0, 00) ; L°°(T)) 
and the conjugation theorem A.l. □ 

Proposition A. 2. 

(i) Let V £ W div (Q) and U £ Wj? (fi). Then UV £ W div (0) and 

div(GVU) = V(G- VU) + U6xv{GV) 

(ii) Let H : WL — > K a Lipschitz function and V £ W d i v {£l). Then 

H(V) £ W div (fl) 

and, almost everywhere 

<Kv(GH(V)) = H'{V)G ■ W + H{V)div(V) 
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Proof, (i) is a consequence of the product of a function in W 1 ' 1 ^, +00) ; i 1 (r)) and a function in 

^•""((O.oo); £°°(r)). 

(ii) Let H and V satisfying the hypothesis. First remark that H being Lipschitz and bounded, 
the function H(V) is in L 1 (fl). Now define V(r,a) = V(<S> T (a)). We will show that H(V)\J<s>\ e 
W^^O.+oo) ; L 1 ^)), in order to apply theorem A.l. The function V is in W^(]0, 00 [ ; L7 oc ( r *)) ( see 
remark 9). Thus it is absolutely continuous in r and _ff being Lipschitz yields H(V) absolutely continuous. 
Hence H{V) E W^QO, 00 [ ; L\ oc {T*)). We conclude the proof by using that 

d T {H{V)\M)=d T {H{V))\M+dYv{G)H{V)\M 

= H\V)d T V\J*\ + &v(G)H(V)\Jg,\ e L^O, oo[xr) 

which requires the following lemma (see [23]) to have d T (H(V)) = H'(V)d T V . 

Lemma A. 2. Let H be a Lipschitz function, I a real interval, X a Banach space and u e W 1 ' 1 ^ ; X). 
Then H ou 6 W 1:1 (L ; X), and almost everywhere 

(Hon)' = H'(u)u' 

□ 
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